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Perturbed version of the complex Toda chain (CTC) has been employed to describe adiabatic inter- 
actions within iV-soliton train of the nonlinear Schrodinger equation (NLS). Perturbations induced 
by weak quadratic and periodic external potentials are studied by both analytical and numerical 
means. It is found that the perturbed CTC adequately models the iV-soliton train dynamics for 
both types of potentials. As an application of the developed theory we consider the dynamics of 
a train of matter - wave solitons confined to a parabolic trap and optical lattice, as well as tilted 
periodic potentials. In the last case we demonstrate that there exist critical values of the strength 
of the linear potential for which one or more localized states can be extracted from a soliton train. 
An analytical expression for these critical strengths for expulsion is also derived. 



I. INTRODUCTION 

One of the most remarkable phenomena occurring in 
nonlinear systems is the possibility to form localized 
states of soliton type as a result of the interplay between 
dispersion and nonlinearity. These states display inter- 
esting properties in response to external fields and their 
interactions have been the subject of continuous inter- 
est since the creation of the soliton theory. Besides the 
motivation from the viewpoint of fundamental physics, 
studies on soliton interactions are very important for ap- 
plications such as optical fiber communications systems, 
where optical solitons are used as information bit carri- 
ers pj. In this context the interaction between neigh- 
boring solitons in a train limits the transmission capac- 
ity of the communication system, so that soliton inter- 
action becomes very important for optimal information 
packing and transmission rates design. Trains of soli- 
tons (fiuxons) in interaction play an important role also 
in long Josephson junctions where their dynamics, in- 
duced by external magnetic fields, is used to construct 
flux-flow oscillators, devices of interest for applications 
in superconducting mm and sub- mm wave electronics . 
Other rapidly developing fields in which interacting soli- 
tons play a crucial role are photonic crystals and Bose- 
Einstein condensates (BEC). 

Recently, BEC solitons have attracted a great deal of 
interest both from the theoretical and the experimen- 
tal point of view (see e.g. reviews |^). In particular, 
self-trapped states capable to propagate in space with- 
out distortion are of interest for pulsed atomic soliton 
lasers, atomic nano-litography and high precision inter- 
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ferometry [3 and the transport of BEC solitons in the 
presence of external potentials, serving as magnetic and 
optical traps and waveguides may become important in 
future technologies. 

The aim of the present paper is to study the adia- 
batic dynamics of a train of N interacting solitons of 
the nonlinear Schrodinger equation (NLS) in weak exter- 
nal potentials. As a physical model we consider matter- 
wave solitons in quasi-one dimensional BEC with attrac- 
tive interactions between atoms, such as 7 Li, 85 Rb or 
137 Cs. The results, however, are of interest also for non- 
linear optics and for photonic crystals. In particular, we 
study the N-soliton dynamics in parabolic, linear and pe- 
riodic potentials, modelling the NLS train soliton inter- 
action in terms of a complex Toda chain (CTC) which 
is valid in the adiabatic approximation. This model 
has been successfully used in previous papers (see Refs. 
[H> \Hi EJ [3 Il9l|23') and will be used here to continue 
the analysis in on the perturbed CTC (PCTC) as 
a model for N-soliton interaction in BEC with quadratic 
and periodic potentials. Our results provide an addi- 
tional confirmation of the stabilization properties of the 
periodic potentials observed in [27ll28| in a different phys- 
ical setup. In particular, for the case of an A-soliton 
train trapped in a weak parabolic trap we find that the 
train performs contracting and expanding oscillations if 
its center of mass coincides with the minimum of the 
potential, while it oscillates around the minimum of the 
potential as a whole if its center of mass is shifted from 
the minimum. In the last case contracting and expanding 
motions of the soliton train is superimposed to the center 
of mass dynamics. As the strength of the parabolic trap 
increases we find from numerical simulations that the N- 
soliton dynamics becomes more complicated with the the 
merging (splitting) of individual solitons when the train 
is contracted (expanded) during its oscillating motion in 
the trap. This behavior resemble the phenomenon of 
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"missing solitons" observed in the experiment gg. 

We have investigated the case of a tilted periodic po- 
tentials i.e. a potential which is the superposition of pe- 
riodic and linear potentials. The effect of the linear po- 
tential, if it is strong enough, is that it can overcome the 
stabilization effect of the periodic potential. As a result 
we show that there exist critical values of the strength 
of the linear potential for which one or more localized 
states can be extracted from a soliton train (array). We 
find that the critical strength for expulsion changes with 
the number of the solitons in the train. To this regard 
we have derives an analytical expression for the poten- 
tial critical strengths at which expulsion is achieved by 
means of an Hamiltonian approach for the PCTC. From 
this comparison we find that the analysis based on the 
PCTC provides a good description for the expulsion phe- 
nomena. 

We remark that since the critical value of the potential 
strength is an indirect measure of the binding energy of 
an N-solitons train, one could use BEC soliton arrays in 
accelerated optical lattices to measure N-solitons matter 
waves binding energies. One could indeed reproduce the 
effect of the linear potential by means of accelerated op- 
tical lattices and measuring the critical acceleration at 
which the expulsion of one soliton from the array occurs. 
We hope experiments in this direction will be performed 
soon. 

The paper is organized as follows. In section II we in- 
troduce the mean field Gross-Pitaevskii equation (GPE) 
appropriate for quasi one dimensional BEC in external 
potentials and discuss the N-soliton interaction in terms 
of the complex Toda chain. In Section III we study the 
perturbed nonlinear Schrodinger equation and the per- 
turbed CTC equation in the adiabatic approximation for 
different types of external potentials: quadratic, linear 
and periodic. The analysis of the N-soliton dynamics 
obtained from the perturbed CTC with the above poten- 
tials is compared with direct numerical GPE simulations 
in section IV. In section V we introduce a Hamiltonian 
formulation of the PCTC and derive an analytical ex- 
pression for the critical strengths of the linear potential 
for which the phenomenon of soliton expulsion occurs. 
Finally, in the last section the main results of the paper 
are briefly summarized. 



II. MODEL EQUATIONS AND N-SOLITON 
INTERACTION 



The dynamics of a condensate in the mean-field ap- 
proximation at zero temperature is governed by the 
3D nonlinear Schrodinger equation (in the BEC context 
called also as Gross-Pitaevskii equation (GPE)) 



- — V 2 + V{x, yi z) 
2m 



Airh a s Af 2 



(2.1) 



where ^(r,^ is the macroscopic wave function of the 
condensate normalized so that J" |\I r (r)| 2 cfr = 1, J\f is the 
total number of atoms, m is the atomic mass, a s is the s- 
wave scattering length (below we shall be concerned with 
attractive BEC for which a s < 0), and 

V(x,y,z) = j[u; 2 x x 2 +u; 2 x (y 2 + z 2 )} (2.2) 

is the axially symmetric trapping potential which pro- 
vides for the tight confinement in the transverse plane 
(y, z), as compared to loose axial trapping, assuming 
uj 2 /oj 2 l <C 1. The condensate trapped in such a potential 
acquires highly elongated form. 

When the transverse confinement is strong enough, so 
that the transverse oscillation quantum, Hu±, is much 
greater than the characteristic mean-field interaction en- 
ergy, (47rfi, 2 a s /m)A/'| x I / | 2 , the dynamics is effectively one- 
dimensional. In this case, the wave function may be ef- 
fectively factorized as ^(x, y, z, t) — ip( x i tyfiiVi z )i where 
4>(y,z) — exp[— (y 2 + z 2 )/2a 2 L )/ y / ira± is the normalized 
ground state of the 2D harmonic oscillator in the trans- 
verse direction, with a± = ^/fij rmJ± being the corre- 
sponding transverse harmonic-oscillator length. Insert- 
ing the factorized expression into the 3D GPE 1|2.1|) , and 
integrating it over the transverse plane (y, z), one derives 
the effective ID equation 



dip 



2m dx 2 



ulx 2 + g 1D N\ijY 



rl>, (2.3) 



where we have neglected the zero-point energy of the 
transverse motion, ?iuj±, and defined a coefficient of the 
ID nonlinearity, gm — ^nh 2 a s m~ Y J \4>(y, z)\ A dydz = 
2a s fiuj±. It is convenient to use normalized units for 
time and space variables, introducing transformations: 
t — ► uj±t, x — > x/ax, and the rescaled wave function 
u ■ 



du 



ld 2 u 
2dx2 



V 2 x 2 u + \u\ 2 u = 0, 



(2.4) 



where V2 — ui 2 / (2uj 2 _) , is a small parameter characteris- 
ing the strength of the external parabolic potential. In 
the experiment of Rice university 6] a matter-wave soli- 
ton train (comprising N ~ 10 solitons) of Bose-condensed 
7 Li atoms (a s = —0.21 nm, m — 11.65 • 10~ 27 kg) was 
created. Radial confinement was strong, uj± ~ 800 Hz, 
while the axial one had been as weak as oj x ~ 3 Hz. In 
the experiment 0, where a single soliton of 7 Li BEC 
was created, the trap aspect ratio was smaller: loj_ ~ 
710 Hz, uj x - 50 Hz. Therefore, V 2 ~ 10" 5 + 10~ 3 is 
in the range of realistic experimental conditions. Below 
we shall consider also other kinds of weak potentials in 
the axial direction x, instead of (or combined with) the 
parabolic trap in Eq. I|2.4II . see 0|2j|. Assuming them 
as perturbations iR[u] = V(x)u(x,t), we move them to 
the right hand side of the governing equation. 

The Af-soliton train interactions for the nonlinear 
Schrodinger equation (NLS) and its perturbed versions 
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.du 1 8 2 U . |2 . 1 fr. r\ 

1 at + 2dx* + |u| M ' ( 5) 

started with the pioneering pap er [ a, b y now has been 
extensively studied (see 0, EU EH Ej> D3 an d refer- 
ences therein). Several other nonlinear evolution equa- 
tions (NLEE) were also studied, among them the modi- 
fied NLS equation El El El El El, some higher NLS 
equations 13], the Ablowitz-Ladik system EH an d oth- 
ers. 

Below we concentrate on the perturbed NLS eq. I|2.5[l . 
By ./V-soliton train we mean a solution of the (perturbed) 
NLS fixed up by the initial condition: 



N 

u(x,t = Q) = Y / u 1 k s {x,t = 0), 

k=l 

z k (x,t) = 2v k (x - &(*)), 
<f) k (x,t) = —z k +S k (t), i 



t 2v k e 1 ^ 

<0M) = — r — . 

cosh z k 

(2.6) 

&(t) = 2A* fc t + &,o (2.7) 
«fc(t)=W fc t + <y fc) o, (2.8) 



Each soliton has four parameters: amplitude v kl velocity 
fi k , center of mass position £ k and phase S k . The adia- 
batic approximation uses as a small parameter Eq <ti 1 
the soliton overlap which falls off exponentially with the 
distance between the solitons. Then the soliton parame- 
ters must satisfy @: 

K--^o|<^o, |m* - Mo I < Mo, 
|»*-»*||&+i,o-&,o|»l, (2.9) 

where v = jfJ2k=i u k> and fi = i X)fc=i arc tlic 
average amplitude and velocity respectively. In fact we 
have two different scales: 

, 1/2 , , 1/2 

H - fol - V i iMfc - Mol - e > 

|Cfe+i,o - £fc,o| — £o ■ 

One can expect that the approximation holds only for 
such times t for which the set of 4iV parameters of the 
soliton train satisfy (|2.9[) . 

Equation 12. 5|) finds a number of applications to non- 
linear optics and for R[u] = is integrable via the in- 
verse scattering transform method |2(J, |2l|. The N- 
soliton train dynamics in the adiabatic approximation is 
modelled by a complex generalization of the Toda chain 

^% = 16^ 2 ( e «i+i-«* - e^-^-A , j = 1, ... ,N. 
d^ 

(2.10) 



The complex- valued Q k are expressed through the soliton 
parameters by: 

Q k {t) = 2i\ t k (t)+2kln(2v )+i(kTr-6 k (t)-5 ), (2.11) 

where 5q = l/JVJ)it=i $k ano - -^o = Mo + w'o- Besides we 
assume free-ends conditions, i.e., e~®° = e® N+1 = 0. 

Note that the /V-soliton train is not an /V-soliton solu- 
tion. The spectral data of the corresponding Lax opera- 
tor L is nontrivial also on the continuous spectrum of L. 
Therefore the analytical results from the soliton theory 
can not be applied. Besides we want to treat solitons 
moving with equal velocities and also to account for the 
effects of possible nonintegrable perturbations R[u]. 

The fact f2 %]x% that the CTC, like the (real) Toda 
chain (RTC) [23], is a completely integrable Hamiltonian 
system allows one to analyze analytically the asymptotic 
behavior of the ./V-soliton trains. However unlike the 
RTC, the CTC has richer variety of dynamical regimes 
HElliil such as: 

• asymptotically free motion if Vj 7^ v k for j 7^ k; this 
is the only dynamical regime possible for RTC; 

• N-s bound state if v\ — ■ ■ ■ — vn but £ k 7^ Q for 

k + j; 

• various intermediate (mixed) regimes; e.g., if V\ = 
V2 > ■ ■ ■ > wjv but Cfe 7^ Cj f° r k 7^ j then we will 
have a bound state of the first two solitons while 
all the others will be asymptotically free; 

• singular and degenerate regimes if two or more of 
the eigenvalues of L become equal, e.g., Ci = C2 • • • 
and (j 7^ ( k for 2 < j 7^ k. 

By Cfc — v k + iwk above we have denoted the eigen- 
values of the Lax matrix L in the Lax representation 
L T = [M, L] of the CTC where: 

N N-l 

L = ^b k E kk + a k (E k , k+1 +E k+hk ), (2.12) 
fe=i fc=i 

where 

1 dQk _ Mfc + iv k _ 1 t Qh+1 -Q h )/2 

h ~ 2 air ~ 2 ' k ~ 2 

and the matrices E kp arc defined by (E kp )ij — S k iS P j. 
The eigenvalues £fc of L are time independent and 
complex- valued along with the first components n k = z± 
of the normalized eigenvectors of L: 

L^=^ k \ (z< k \z< m ^ = S km . (2.13) 

The set of {£ k = v k + iw kl rj k = a k +i9 k } may be viewed 
as the set of action-angle variables of the CTC. 

Using the CTC model one can determine the asymp- 
totic regime of the iV-soliton train. Given the initial pa- 
rameters M/b(0) 5 ^(0), £fc(0)i S k (Q) of the iV-soliton train 
one can calculate the matrix elements b k and a k of L 
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at t = 0. Then solving the characteristic equation on 
L\t=o one can calculate the eigenvalues Oc to determine 
the asymptotic regime of the ./V-soliton train 0, • An- 
other option is to impose on a specific constraint, e.g. 
that all £fc be purely imaginary, i.e. all Vk = 0. This will 
provide a set of algebraic conditions L\t=o, and on the 
initial soliton parameters /U/s(Q), ^fe(O), £fc(0), 5k(Q) which 
characterize the region in the soliton parameter space re- 
sponsible for the iV-soliton bound states. 



III. THE PERTURBED NLS AND PERTURBED 
CTC EQUATIONS 

We will consider several specific choices [u] of per- 
turbations, p = 1,2, . . . in H2.5fl . In the adiabatic approx- 
imation the dynamics of the soliton parameters can be 
determined by the system (see @ for N = 2 and 0, 
for N > 2): 



~dT 
dt 



-4z/ (e Wh 



Qk+i-Qk 



a Qk-Qk- 







(p) 



iN, 



(p) 



(3.1) 



d&k 
~~dt 



2(^ + ^)+4%2) 



where Afc = /ife + ii/fc and X 



(v) 



On WW 



D 



(p) 



The 



right hand sides of Eqs. I|3.1[l - H3.2[) are determined by 
ftjf ) [u] through: 



(p) 



N, 



~(p) 



1 



2 7_oo coshz/c 
1 f°° dzk sinh Zfc 
2 



^Re (4 P) Me- i0 * 



1 



cosh Zfe 



Im JjW[«] ( 



(3.3) 
,(3.4) 



4ffc 7-oo coshzfc 



Re (i#>[ u ] e -**») , (3.5) 



(p) 



1 



(izfe (1 — Zfc tanh z^) 



Im (i?i p) Ne- 



2 ^fc 7-oo cosh^ 
Inserting QH.lj) . I|3.2(l into (|2.11|) we derive: 



(3.6) 



(P) y{p) ^(p^ 



iV 



(3.7) 



N L^i ■ 



AT ^ J 
i=i 



(p) 



3=1 



In deriving eq. I|3.7|l we have kept terms of the order 
Ai/k ~ ©(y/eo) and neglected terms of the order O(eo). 



The perturbations result in that vq and /^o may become 
time-dependent. Indeed, from (|3.1|) we get: 



^0 _ . ,(p) 



(3.8) 



The small parameter eo can be related to the initial 
distance tq = \& — £,i\t=o between the two solitons. As- 
suming v\ t 2 ~ j/q we find: 



eo 



dx \u\ s {x, Q)u^{x, 0) I ~ 8^ r e 



-2u r Q 



(3.9) 



In particular, (|3.9|l means that eo ~ 0.01 for ro ~ 8 and 
z/q = 1/2. 

We assume that initially the solitons are ordered in 
such a way that £fc+i — ~ ro- One can check [ToL ITflj 



that N, 



(p) 



M, 



(p) 



exp(— 2vo\k — p|ro). Therefore the 
interaction terms between the fc-th and fc ± 1-st solitons 
will be of the order of e - 2v ° r a^ the interactions between 
fc-th and k ± 2-nd soliton will of the order of e~ Al, ° r ° <C 

e -2u r ^ 

The terms Sj^ , are of the order of 

exp(— 2^o ro), where a = or 1. However they 
can be neglected as compared to Jlk and where 



eo, (3.10) 



(p) 

The corrections to Nl , . . . , coming from the terms 
linear in u depend only on the parameters of the fc-th 
soliton; i.e., they are 'local' in fc. The nonlinear terms in 
u present in iR^ [u] produce also 'non-local' in fc terms 



in AT, 



(p) 

k ' 



A. Quadratic and tilted potentials 

Let iR[u] — V(x)u(x,t). Our first choice for V(x) is a 
quadratic one: 



V {1) {x) = V 2 x 2 + V lX + Vb- 
Skipping the details we get the results: 



(i) _ 

k " u ' 



M, 



(i) 



57(1) - 

-^2& ~ y, 



(3.11) 

(3.12a) 
(3.12b) 



W = V2 



2 -€fc)-Vi^-Vo, (3.12c) 



and — As a result the corresponding PCTC 

takes the form Hal: 



dXk 
dt 



V 2 tk - y, 



(3.13) 



dQ k 
dt 



. N 



3=1 
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where \ k —^k + i^k- 

If we now differentiate (|3.14[1 and make use of (|3.13|) 
we get |T^ : 



^ = 16, 2 (e^- 
+ 4^ ( V 2 £ k + 



Qk—Qk-i 



(3.15) 



Vi 



dt N ^— ' dt 



It is reasonable to assume that V 2 ~ O(eo/N); this 
ensures the possibility to have the A-soliton train 'inside' 
the potential. It also means that both the exponential 
terms and the correction terms are of the same 

order of magnitude. From eqs. i|3.13|l and (|3.14|1 there 
follows that dvQ /dt = and: 



"^--^o-y, -^ 



(3.16) 



where /io is the average velocity and £,o — jj £j) ^ s 

the center of mass of the ./V-soliton train. The system of 
equations (|3.16(l for V 2 > has a simple solution 



/■to(f) = Moocos($(t)), 
~2~ 



-/zoo sin($(t)) 



A. 

2F 2 : 



(3.17) 



where $(t) = v / 2T^i+ < I>o, and //oo and <&o are constants of 
integration. Therefore the overall effect of such quadratic 
potential will be to induce a slow periodic motion of the 
train as a whole. 

By tilted potential below we mean a particular case of 
the quadratic potential with V 2 — 0. Then the equations 
(j3~T^I - (j3~T7|l take the form: 



d^Q Vi d£ _ ~ 

and have the simple solution 

Mow = — jt + ii Q o, 



Hoot + Coo, 



(3.19) 



Therefore the tilted potential accelerates the soliton train 
as a whole in a prescribed direction. It can be used also 
to 'pick up' and accelerate one or more of the solitons 
from the train confined to a periodic potential. 



B. Periodic potentials 

One may consider several physically important choices 
of periodic potentials. The simplest one is 

V (2) (x) = A x cos(Oia; + Q ) 

= Ax- 2Ax sxn 2 (ttxx/2 + fi /2), (3.20) 



where A, f2 and f^o are appropriately chosen constants. 
NLS equation with similar potentials appear in a natural 
way in the study of Bose-Einstein condensates, see [24J . 

The relevant integrals for N k , Mk, Sfe and Dk are equal 
to El: 



N, 



(2) 



"(2) _ 
^k — u > 



M, 



(2) 



0, 

nAxflf 1 



Dl 2) = - 



8v k sinh Z k 
tt 2 Ax^I coshZ fe 



sin(Oi^ fc + Q ), 



(3.21) 
(3.22) 



■ 7 C0S ("A + fio), (3-23) 
to^fe sinn Zk 



where Zx,k = 7i"^i/(4^fc). These results allow one to de- 
rive the corresponding perturbed CTC models. Again we 
find that dv^/dt = 0. 

The second case we will consider is a linear combina- 
tion of two potentials of the form (|3.20() with correlated 
frequences: 

V {3) (x) = Ax cos(fiia; + fi ) + M cos(n 2 x + Q ) (3.24) 

where Aj, Qj and f2o are appropriately chosen con- 
stants. Such potentials with correlated frequences, e.g. 

= 2f2i also _appear in the study of Bose-Einstein con- 
densates, see [24). 

The relevant integrals for Nk, Mk, and Dk are equal 
to El: 



ivf = 0, 
r(3) 



- 



(3.25) 



Ml ' = AxMk(nx,Zx,k)sm(rLxZk + n a ) 



+ A 2 M k (n 2 ,Z 2 ,k)sin(n 2 Z k +n ), (3.26) 

= AxD k (ni,z lik )cos(n^ k + n ) 

+ A 2 D k (n 2 ,z 2 , k )cos(n 2 z k + n ), (3.27) 



(3.18) where 



M k {Q, h Z jtk ) 
D k (Clj, Zj t k) 



1 



lv k sinh^fe' 
tt 2 ^ 2 coshes 



16^ S inh 2 Z 



(3.28) 



Zj.k = nflj/(4vk)- These results allow one to derive the 
corresponding perturbed CTC models. Again we find 
that dvo/dt = 0. 

The last case we consider is an elliptic potential of the 
form: 



V (4) (af) = -Bsn 2 (nx;k) 

OO 

= -B2~2G s (k) sin{n s x), 



G s (k) = 



(2S + 1) 3 7T 2 



2fc 3 



8fc 3 K 2 J Ksmh(w a ] 
(2s + 1)ttK' 



(3.29) 



0.30) 



2K 



(3.31) 
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where k is the module of the elliptic function, K = K(k) 
is the complete elliptic integral of the first kind, and B 
is a constant. 

The second line of formula H3.29[l provides the expan- 
sion of sn 2 (fix; k) as infinite series of trigonometric func- 
tions. Each of the terms in this series can be treated as 
perturbation just like above assuming Qq = —tt/2. As a 
result we get: 



N, 



(4) 



0. 



"(4) 



0. 



(3.32) 



The evolution of a -/V-soliton train in the absence of po- 
tential (V(x) = 0) is well known, see e.g. 0,0. These 
papers propose a method to determine the asymptotic 
dynamical regime of the CTC for a given set of initial 
parameters ^fc(O), ^fc(O), £fc(0) and 4(0). Below we will 
use mainly the following set of parameters ^fc(O) = 1/2, 
/Ufe(0) = 0, £fc+i(0) — £fc(0) = r with two different choices 
for the phases: 



4(0) = kn, 

4(0) = o. 



(4.1) 
(4.2) 



M< 4) = B^M k (n s ,Z S!k )G s {k)cos(n s Z k ), (3.33) 



s=0 



D k = -Bj2 D k(n s ,Z Stk )G s {k)sm({l s Z k ), (3.34) 



s=0 



where M k , D k and f2 s = (2s + l)f2 are defined as 



IV. CTC ANALYSIS AND COMPARISON WITH 
NUMERICAL SIMULATIONS 

The dynamics of an individual soliton in a train is de- 
termined by the combined action of external potential 
and the influence of neighboring solitons. The interac- 
tion with neighboring solitons can be either repulsive, 
or attractive depending on the phase relations between 
them. Particularly, if their amplitudes are equal and the 
initial phase difference between neighboring solitons is 7r 
(as considered below) they repel each other giving rise 
to expanding motion in the absence of an external field 

urn 

The external potential counterbalances the expansion, 
trying to confine solitons in the minima of the potential. 
It is the interplay of these two factors - the interaction of 
solitons and the action of the external potential, which 
gives rise to a rich dynamics of the ./V-soliton train. 

To verify the adequacy of the perturbed CTC model 
for the description of the ./V-soliton train dynamics in 
external potentials we performed comparison of predic- 
tions of corresponding perturbed CTC (PCTC) system 
and direct simulations of the underlying NLS equation 
(|2.5|) . Below we present results pertaining to a matter- 
wave soliton train in a confining (i) parabolic trap and in 
(ii) a periodic potential modelling an optical lattice. 

Here we present the numerical verification of the 
PCTC model. The perturbed NLS eq. is solved 

by the operator splitting procedure using the fast Fourier 
transform puj . In the course of time evolution we mon- 
itor the conservation of the norm and energy of the N- 
soliton train. The corresponding PCTC equations are 
solved by the Runge-Kutta scheme with the adaptive 
stepsize control |32T |. 



These two types of initial conditions (IC) are most widely 
used in numeric simulations. 

In the absence of potential the IC <|4.1[l ensure the so 
called free asymptotic regime, i.e. each soliton develops 
its own velocity and the distance between the neighbor- 
ing solitons increases linearly in time. At the same time 
the center of mass of the soliton train stays at rest (see 
the left panel of Fig. [TJ. Under the IC (|4.2|l the solitons 
attract each other going into collisions whenever the dis- 
tance between them is not large enough. 

From mathematical point of view the IC l|4.1l) reduce 
the CTC into a standard (real) Toda chain for which the 
free asymptotic regime is the only possible asymptotical 
regime. On the contrary the IC l|4.2|) lead to singular 
solutions for the CTC (see H El Hfl). The singulari- 
ties of the exact solutions for the CTC coincide with the 
positions of the collisions. 

Below we will study the effects of the potentials for 
both types of IC. One may expect that the quadratic 
potential will prevent free asymptotic regime of IC 14.1(1 
no matter how small V2 > is and would not be able to 
prevent collisions in the case of IC 14. 2|) . The periodic 
potential, if strong enough, should be able to stabilize 
and bring to bound states both types of IC. 



A. Quadratic potential 

For the quadratic external potential V(x) = V2X 2 + 
Vix + Vo the perturbed CTC equations in terms of soliton 
parameters have the form: 



dfJ, k 
dt 



dv k 
dt 



e -2Mi*-tk-x) GOS^i) - V 2 tk - y, (4.3) 



!tt = 2 ^ 



dt u k 



; -2^ (? fc -a-i) sin $ fc _ 1 j ; (4.4) 

(4.5) 



V 2 



48za 2 



8 )-V 1 ^ k -Vo, (4.6) 



7 



- 'M - VH l M'- 1 ^ 60: 




jULM-jM-AXJL 
JULiJL^juUJL 

jlaj_«jkaaaj 

_LL*JU4jULS 

iiJJUAiJul 



A k A. A | A A A A 



!mmma!A_!_ 

; ; ' GC ' _ 
ffiuJ '*■)'■■ '■■ 



^ A A A A /\ A A i 






l _* 

: 



r At'* 



-8D -60 -40 -20 



20 40 60 SO 




FIG. 1: Left panel: 9-soliton train with initial parameters as 
in eq. 114. in with ro = 8 in the absence of a potential goes 
into free asymptotic regime. Solid lines: direct numerical 
simulation of the NLS equation 12.51 : dashed lines: as 
predicted by the CTC equations ijl^ - ljUTt with Vo = Vi = 
V2 = and ro = 8. Right panel: Evolution of a 9-soliton train 
with the same initial parameters in the quadratic potential 
V(x) = V2X 2 with V2 = 0.00005. Solid lines: direct numerical 
simulation of the NLS equation 12.51 : dashed lines: solution of 
the PCTC equations 114. 31 - 114.61 . Initially the train is placed 
symmetrically relative to the minimum of the potential at 
x = 0. 



$ fc = 2jU (6»+i - 60 + $k - 



(4.7) 



where ^k, ^fc, £/c and 5k for fc = 1, ...,N are the AN 
soliton parameters, see eqs. H2.6IH2.7fl . 

The effect of the quadratic potential on the ./V-soliton 
train with parameters (|4.1|l is to balance the repulsive 
interaction between the solitons, so that they remain 
bounded by the potential, as illustrated in the right pan- 
els of the figures and [21 The quadratic potentials are 
supposed to be weak, i.e. we choose V% so that 



(4.8) 



Figures and [21 show good agreement between the 
PCTC model and the numerical solution of the perturbed 
NLS equation (|2.5() . They also show two types of effects 
of the quadratic potential on the motion of the 7V-soliton 
train: (i) the train performs contracting and expanding 
oscillations if its center of mass coincides with the mini- 
mum of the potential, (ii) the train oscillates around the 
minimum of the potential as a whole if its center of mass 
is shifted. In the last case contracting and expanding mo- 
tions of the soliton train is superimposed to the center 
of mass dynamics. As one can see from the figures the 
period of this motion matches very well the one predicted 
by formula 1)3.17)1 . Indeed, from eq. (|3.17f) it follows that 
the period of the center of mass motion is T — 2tt/\/2\^. 
For the parameters in fig. [21 we have T ~ 628 (for 9- 
soliton train), T ~ 140 (for 3-soliton train). Similar is 
the dynamics also for the 7-soliton train on Fig. [21 for 
the parameters choosen there we have T ~ 314, in good 



FIG. 2: Harmonic oscillations of a N-soliton train initially 
shifted relative to the minimum of the quadratic potential 
V(x) = V 2 x 2 . Left panel: 9-soliton train, V 2 = 0.00005. 
Right panel: 3-soliton train, V 2 = 0.001. The IC of the both 
trains is given by 1)4. l|l with ro = 8. In both panels solid lines 
correspond to direct simulations of the NLS equation 12,51 . 
and dashed lines to numerical solution of the PCTC equations 
llOll - Oil . 
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FIG. 3: Dynamics of a 7-soliton train placed asymmetrically 
relative to the minimum of the trap V{x) = 0.0002j: 2 . Solid 
lines: results of direct numerical simulations of the NLS equa- 
tion 12.51 . Dashed lines: result of solution of the PCTC sys- 
tem 14.31 - 14.61 for the center of mass f j . The parameters of 
solitons are the same as in 1)4. 1^ with ro = 8. The initial shift 
of the soliton train relative to the minimum of the parabolic 
trap is 10-7T. 



agreement with the numerical simulations. The direct 
simulations of the NLS equation 1)2.5)1 shows that stronger 
parabolic trap may cause merging of individual solitons 
at times of contraction, and restoring of the original con- 
figuration when the train is expanded. This behavior 
reminds the phenomenon of " missing solitons" observed 
in the experiment |6j. However, this situation is beyond 
the validity of the PCTC approach. 
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B. Tilted periodic potential 

Now we consider the dynamics of a N - soliton train in 
a tilted periodic potential, which is the combination of 
periodic and linear potentials 



V(x) = Acos(ilx + too) + Bx. 



(4.9) 



This potential is of particular interest in studies of Bosc- 
Einstein condensates. A train of repulsive BEC loaded 
in such a potential (where the periodic potential was a 
ID optical lattice and the linear one was due to the grav- 
itation) exhibited Bloch oscillations • At each period 
of these oscillations condensate atoms residing in indi- 
vidual optical lattice cells coherently tunneled through 
the potential barriers. This was the first experimental 
demonstration of a pulsed atomic laser 33]. Recently 
a new model of a pulsed atomic laser was theoretically 
developed in (34|, where the solitons of attractive BEC 
were considered as carriers of coherent atomic pulses. 

Controlled manipulation with matter - wave solitons is 
important issue in these applications. Below we demon- 
strate that solitons of attractive BEC confined in optical 
lattice can be flexible manipulated by adjustment of the 
strength of the linear potential. In Fig. 0|we show the ex- 
traction of different number of solitons from the 5-soliton 
train by increasing the strength of the linear potential B, 
as obtained from direct simulations of the NLS equation 
(1) and numerical integration of the PCTC system (38) 
- (41) with 



D 



(2) _ 
k ~ 

(2) _ 



irAn 2 



1 



1 



. sm(n$ k + n )--B,(4.10) 
sinn Zk 2 



it 2 AO? coshZ fe 



16i/| sinh 2 Z fe 



cos(f2£fc + f2o) 



(4.11) 



As is evident from Fig. 0] the PCTC model provides 
adequate description of the dynamics of a N-soliton train 
in a tilted periodic potential. A small divergence between 
predictions for the trajectory of the left border soliton in 
Fig. 0| (d) , is due to the imperfect absorption of solitons 
from the right end of the integration domain. Reflected 
waves enter the integration domain and interact with soli- 
tons, which causes the discrepancy. 



C. Periodic potential 

Another external potential in which the A^-soliton train 
exhibits interesting dynamics is the periodic potential of 
the form V{x) = Acos(Vlx + £Iq). This case also may 
have a direct relevance to matter - wave soliton trains 
confined to optical lattices. The PCTC system in terms 
of soliton parameters has the form: 
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FIG. 4: Controlled withdrawal of solitons from the 5-soliton 
train by adjusting the strength of the linear potential (14.91 
with parameters: A = -0.0005, Q = 2tt/9, H = 0. De- 
pending on the tilt, different number of solitons can be pulled 
out of the train: (a) one soliton at B — —0.00003, (b) two 
at B = -0.00011, (c) four at B = -0.0002 and (d) five at 
B — —0.0003. The initial phase difference and separation 
between neighboring solitons in the train are, respectively, n 
and 9. Initially the train is shifted by — 107r with respect to 
x — for graphical convenience. Solid and dashed lines corre- 
spond, respectively, to direct simulations of the NLS equation 
12.51 and numerical integration of the PCTC system 14.121 - 



dvk 
~~dt 



dt 
dS k 

~dt 



e -2n,K*-e*-0 cos $ fe _i) + M< 2) K), (4.12) 
r 2 M<Lh+i-£k) sm $ fc 

e -2 I ,o&-? fc -i) sin$fc _ i 



(4.13) 
(4.14) 
(4.15) 



where M^\v k ), D h 2) [ Vk) are given in l|3.21|l and i|3.23[l . 
and $fc - in eq. [|4.7fl . 

Each soliton of the train experience confining force of 
the periodic potential and repulsive force of neighboring 
solitons. Therefore, equilibrium positions of solitons do 
not coincide with the minima of the periodic potential. 
Solitons placed initially at minima of the periodic poten- 
tial (Fig. |SJ) perform small amplitude oscillations around 
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FIG. 5: Solitons (continuous line) remain confined around the 
minima of the periodic potential V(x) — Acos(x) (dashed 
line) performing small amplitude oscillations if its strength is 
big enough A — —0.1. 
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FIG. 6: Left panel: The expulsion of solitons from the train, 
as obtained from direct simulations of the NLS equation 12.5|l 
(solid lines), and as predicted by PCTC system H4.121 - H4.15ll 
for the center of mass £i (dashed lines). The IC of the 7- 
soliton train are given by 14. H with ro = 8; the parameters 
of the periodic potential V(x) = Acos(Q.x + f2o) are A = 
—0.001, Q = 7r/4, ilo = 0. Right panel: Oscillations of 
the 5-soliton train with IC given by 14.11 with ro = 9; in a 
moderately weak periodic potential, A — —0.0005, Q = 2tv/9, 
Qo — 0. Solid and dashed lines correspond, respectively, to 
numerical solution of the NLS eq. 12.51 and PCTC system 

ism - 



potential wells under repulsive forces from neighboring 
solitons. When the amplitude of oscillations of particu- 
lar solitons grow and two solitons closely approach each 
other, a strong recoil momentum can cause the soliton to 
leave the train, overcoming barriers of the periodic po- 
tential. In Fig. this happens with bordering solitons 
(the other solitons remain bounded under long time evo- 
lution). It is noteworthy to stress that this phenomenon 
is well described by the PCTC model, as is evident from 
Fig. IHI left panel. 

On the right panel of the same figure we have similar IC 
as in (|4.1I) and we have choosen again the initial positions 
of the solitons to coincide with the minima of the periodic 
potential V(x) — Acos(£lx + fio); i.e. ro = 27r/f2. The 
values of A = —0.0005 and ro = 9 in the right panel 
of Fig. [f)| now are such that the solitons form a bound 
state. Therefore for any given initial distance ro there is 
a critical value A cr (ro) for A such that for A > A cr (ro) 
the soliton train with IC l|4.1|l will form a bound state. 

In contrast to the quadratic potentials, the weak pe- 
riodic potential is unable to confine solitons, and re- 
pulsive forces between neighboring solitons (at z/jfc(0) = 
1/2, <5fc(0) = kit) induces unbounded expansion of the 
train similar to what was shown in the left panel of Fig. 
□ 

The periodic potential can play stabilizing role also 
for the IC (|4.2|l . when the zero phase difference between 
neighboring solitons correspond to their mutual attrac- 
tion. If the periodic potential is strong enough, solitons 
do not experience collision. The weak periodic potential 
cannot prevent solitons from collisions, which eventually 
leads to destruction of the soliton train, as illustrated in 
Fig. [7] Again for any given initial distance ro there will 
be a critical value A' CT (r ) for A such that for A > A^. r (ro) 
the soliton train with IC 14. 2|) will form a bound state 
avoiding collisions. 

Attractive interactions at zero phase difference be- 
tween neighboring solitons can be balanced by expulsive 
force on solitons, if the train is positioned on an inverted 
parabolic trap. In this case solitons far from the center 
experience stronger expulsive force and leave the train, 
as illustrated in Fig. 



these minima, provided that the strength of the potential 
is big enough to keep solitons confined. As opposed, the 
weak periodic potential is unable to confine solitons, and 
repulsive forces between neighboring solitons (at phase 
difference it) induces unbounded expansion of the train. 
In the intermediate region, when the confining force of 
the periodic potential is comparable with the repulsive 
forces of neighboring solitons, interesting dynamics can 
be observed such as the expulsion of bordering solitons 
from the train, as shown in the left panel of Fig. [5] This 
phenomenon, revealing the complexity of the internal dy- 
namics of the train, can be explained as follows. Each 
soliton performs nonlinear oscillations within individual 



V. HAMILTONIAN APPROACH TO 
PERTURBED NLS AND CTC 

The Hamiltonian method has played important role in 
the analysis of integrable and close to integrable nonlin- 
ear evolution equations and dynamical systems, see [2l| . 
In this section we will outline how this method can be 
used for the analysis of the A^-soliton interactions. 

The relation between the Hamiltonian properties of the 
NLS eq. and the CTC model was derived in (35|. Here 
we will show that this approach can be extended also to 
the perturbed versions of NLS and CTC. Indeed, the 
Hamiltonian of eq. (|2.5|l with iR[u] = V(x)u(x, t) is 



10 



A, 



A A A A A 

A A A A. A 
A AAA A. 
A A 4 iA 

_A_A A A a 

_AA-A-AA 

AAA A A 
A A A A. A 
AAA A. A 

. A 




-40 



A AAA A 
A A A A A 
AA A A A 
.A a A A , A 



600 



vlAV 



_7'vA...A_ a\_ 
AA A AA 



_a_AAA_ju 

_A_iUA-A_ 

_a A a A_a_. 
_ A A A 55 a 
js\A A A a 

AAJL_AA_ 

1AJUMl_ 

_AJuULA_ 
_JuLAAA 
_A A A A A— 
AA1AA 

A f\ -\ /!'■ 



■150 



ZOO 



100 
t 



■50 



-20 
X 



20 40 



FIG. 7: Dynamics of a 5-soliton train with zero phase differ- 
ence between neighboring solitons, in the periodic potential 
V(x) = Acos^a; + f2 ) with f2 = 7r/4, Q = 0, and r = 8. 
Left panel: When the periodic potential is strong enough 
A = —0.02, the N-soliton train remain confined, each soliton 
performing small amplitude oscillations around the minima 
of individual cells. Right panel: Weaker periodic potential 
A = —0.01 cannot prevent solitons from collisions, which de- 
stroy the train. Solid and dashed lines correspond, respec- 
tively, to numerical solution of the NLS eq. (12.511 and PCTC 
system lETHl - lETTSTl . 
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FIG. 8: Stabilization of a soliton train in a combined poten- 
tial (periodic + inverted parabola): V(x) = — 0.02 cos(f2x) — 
0.000145a; 2 . Separation between in-phase (Sk = 0) solitons is 
ro = 9, Q — 27r/9. Bordering solitons leave the train as they 
experience stronger expulsion, while the central ones remain 
bounded. 



equal to: 



H — i?NLS + Hy, 



H v = 



poo -, 

/ dx-(\u x \ 2 -\u(x,t)\ 4 
dxV{x)\u(x,t)\ 2 , 



(5.1) 
(5.2) 

(5.3) 



One of the ways to derive the CTC from the NLS 
is based on the use of the variational method [TTl l36| . 
Namely one constructs the Lagrangian of NLS, then takes 



an anzatz of the form: 



N 



(5.4) 



fe=i 



and integrate over x neglecting terms of order e k with 
k > 1. Here uj?~ (x,t) is the one-soliton pulse with pa- 
rameters /ife, i>k, £,k and Sk- The results must depend only 
on the AN soliton parameters. It was pointed out in [3f| 
that if we apply this method directly to the Hamiltonian 
A^nls w e get additional singular in e elements which are 
taken care of by a proper regularization. The regularized 
Hamiltonian is 



H Ieg = A{^1 + t$)d - 4/j C* 2 + Hnls, 

poo 

Ci = dx\u(x,t)\ 2 , 



(5.5) 



C 2 



dx —(u x u(x, t) — u*(x, t)u x ), 



Then one can show that 

-ffrcg = 32/yo-ffcTC + const 



(5.6) 



where Hqtc is the Hamiltonian of the CTC. It is ob- 
tained from the Toda chain Hamiltonian: 



N JV-1 



(5.7) 



k=l 



k=l 



by a complexification procedure after which the dynam- 
ical variables become complex- valued: 

Pk -> Pk= Po,k+ipi,k, qk -> Qk = qo,k+iqi,k, (5.8) 
which must satisfy the Poisson brackets: 



{po,k,qo, s } = 4 S , 
{Po,k,qi,s} 



{pi,k,qi, s } = S ks , (5.9) 
0, {/'i i • .- ! (5.10) 



Then the CTC can be written down as a standard Hamil- 
tonian system with 2iV-degrees of freedom and Hamilto- 
nian provided by the real part of the complexified ifrc : 



N 1 

-ffcTC = ?(P0,k-Pl,k) 



fc=l 



JV-1 



16^ J2 e qo - k+1 ~ qo - k cos( 9 i, fc+ i - gi,fc!j5.11) 



fc=i 



Po,k 



dqo,, 



dt 



Pi,k 



dq 



i,k 



dt 



-Av v k . 
(5.12) 

It remains to replace go.fc and q\^ in terms of the soliton 
parameters in order to get the final expression for -ffcTC ■ 
Let us now derive the Hamiltonian for the perturbed 
CTC. To this end we will evaluate Hy in terms of the 
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soliton parameters. Inserting the anzats l|5.4|l into the 
integrand for Hy we obtain two types of terms. The first 
type is 

JV JV 

k,k+l), (5.13) 

fe=l fc=l 

/oo 
dxV{x)\u[ ls \x,t)\ 2 , (5.14) 
-OO 
POO 



In what follows we will neglect the terms Hk,k—i a s 
compared to because their ratio is of the order of 
e. The integrals in l|5.14|l for the quadratic and periodic 
potentials have the form: 



ttAQ 



481/2 



V^k + V 



sinh Z 



k 



cos(f2£ fc + fio)< 



(5.15) 



One can also evaluate the Poisson brackets between the 
soliton parameters inserting the expressions for p a ^ and 
qp tS with a, f3 = 0, 1 into l|5.9|l . We skip the lengthy cal- 
culations here and only note that these Poisson brackets 
combined with the Hamiltonian 



JV 

-Hpctc = Hqtc + 2^ ^ k 
fc=l 

= Hcm + flVj 



(5.16) 
(5.17) 



indeed produce the equations of motion for the PCTC. 

Our next step is to separate the center of mass motion 
described by Hcm- 



Hcm = N (8^1 - i%) + 4i^(V 2 €o + Vi£o + Vb) 



n cos(fi£ + ^o) 



(5.18) 



4i/ ' 



where the subscript stands for the average value of the 
corresponding parameter, e.g. £o = l/A^X^feLi Then 
the Hamiltonian iTy would describe the relative motion 
of the solitons. We will express it as a function of the 
averaged parameters fio, . . . , Sq and the relative param- 
eters: 

Mfc = Mfe - Mo, v k = v k - u , 

« = (5-19) 

Note, that only jV — 1 of the relative parameters are inde- 
pendent; obviously they satisfy 2feLi -^fc = where Xfc 
stands for jufe, . . . , 8 k - 

In evaluating _ffy we will neglect higher order terms, 
i.e. terms of the order e 3 / 2 and higher, as well as terms 



of the order V^e 1 / 2 , Ae 1 / 2 and higher. As a result Hqm 
and i?v simplify to: 



H CM = 8vZ(j4 - i%) + 4z/ I + +V ' ' 



V 2 n 2 



48^ 



nAft 



cos(fi£ + fio), 



sinh Zo 



(5.20) 



fe=i 

JV-l 



,<70,fc+l — <?0,fc 



cos(5i,fe+i - qi,k), (5.21) 



fc=i 



= 4n,V 2 (& - Co) 2 

+ . - - (cos(fi& + Oo) - cos(^ + n )). (5.22) 
smn Zo 

The Hamiltonian i?cM which describes the center of 
mass motion is more simple than Hy and often one is 
able to solve explicitly the corresponding equations of 
motion, see Section IV. The next step would be, using 
the known expressions for the averaged variables /xo(i), 
. . . , 8o(t) to insert them into Hy and try to analyze the 
corresponding equations of motion. This would provide 
us with information about the relative motion of the soli- 
tons around the center of mass. Usually we get a set of 
nonlinear and non-integrable ODE. 

Our idea here is to use the explicit form of Hy for the 
estimation of the critical values of potential strengths A, 
V2, V\ for which the soliton motion becomes qualitatively 
different. Doing this we are making use of the follow- 
ing hypothesis based on the well known fact that bound 
states have negative energies while asymptotically free 
motions should correspond to positive energies. There- 
fore we evaluate the Hamiltonian Hy inserting in it the 
initial soliton parameters along with the known expres- 
sions for /Lto, ■ ■ • , 60 . If the result is negative we may 
expect that the relative motion of the solitons will be 
bounded; otherwise one may expect that at least one (or 
more) of the solitons will move away from the others. 
The critical value of the corresponding constants will be 
derived below with the condition Hy = 0. 

Assume we have only periodic potential present and 
the initial soliton configuration is l|4. lfl . For A = the 
solitons will go into asymptotically free regime. Switch- 
ing on the self-consistent periodic potentials (such that 
the solitons initially are located at its minima) it is nat- 
ural to expect that for A > A cl - the solitons will be stabi- 
lized into a bound state. Then from the condition Hy = 
we get: 



A r 



1 

N 



6 M e -2uoro shlh ^ (5 . 23) 

iril 4^o 



Note that the critical values generically should depend 
not only on the number of solitons N, but also on the 
initial configuration. The approach we used is not very 
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jcxp ; At 



2tt -0.0053 -0.00365 1.45 

7 -0.0025 -0.00166 1.51 

8 -0.00084 -0.00057 1.47 

9 -0.00030 -0.00020 1.50 

TABLE I: The values of AJ: xp obtained from numeric simu- 



lations with PCTC versus A. 
N = 3. 



th 



obtained from eq. I|5.23|l for 



potentials). Both the analytical treatment in the frame- 
work of the PCTC model, and numerical analysis by di- 
rect simulations of the underlying NLS equation show 
that the PCTC is adequate for description of the adia- 
batic JV-soliton interactions in weak external potentials. 
Hamiltonian approach for the perturbed NLS and CTC 
has been developed, and applied to the analysis of the 
soliton "expulsion" from the train, which is confined to 
a periodic potential. 



ro 


A 5,exp 


a 3,cxp 


4 th 


xcxp ; ith 




2tt 


-0.0068 


-0.0029 


-0.0044 


1.55 


0.66 


7 


-0.0031 


-0.0013 


-0.0020 


1.56 


0.65 


8 


-0.00115 


-0.00043 


-0.00068 


1.68 


0.63 


9 


-0.00041 


-0.000155 


- 0.00024 


1.71 


0.65 



MAX 



TABLE II: The values of A^ cxp and Acf Kp obtained from 
numeric simulations with PCTC versus A*5? obtained from 
eq. KZSl for TV = 5. 



sensitive to this. It can not provide us with the interme- 
diary critical values when the soliton train is stabilized 
after emmiting two or more solitons. 

We compared the theoretical predictions for A CI from 
eq. (|5.23|) with the data coming from the numeric solu- 
tions of the corresponding PCTC for different choices of 
the initial distance between the solitons ro. The results 
are collected into the table [I] The conclusion is that eq. 
(|5.23l) provides correct dependence of A cr on ro up to an 
overall constant factor of the order of 1.5. 

In the next table [H] we summarize two sets of critical 
values for the 5-soliton trains. From our numeric exper- 
iment we derive two critical values of A. The first one 
A^ exp describes the value of A above which all 5 solitons 
form a bound state; the second one ^4^. oxp shows the value 
of A above which the three middle solitons form a bound 
state while the two end ones separate. 

Again we see, that formula (|5.23l) describes correctly 
the dependence of both critical values on ro. 



VI. CONCLUSIONS 

We have studied the dynamics of a TV-soliton train 
confined to external fields (quadratic, periodic and tilted 



As a physical system of direct relevance we have consid- 
ered matter- wave soliton trains in magnetic traps and op- 
tical lattices. In the range of parameters for the trapping 
potential, used in the BEC soliton train experiments, we 
found a good agreement between the analytical estimates 
based on the PCTC model and numerical simulations of 
the governing NLS equation. In what concerns the crit- 
ical strength of the periodic potential at which the soli- 
ton expulsion occurs, analytical predictions qualitatively 
agree with numerical simulations. The developed theory 
can be useful for controlled manipulation with matter- 
wave soliton trains. 
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